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THIRD ORDER EQUIVALENT EQUATION FOR THE RELATIVE VELOCITY 
LATTICE BOLTZMANN SCHEMES WITH ONE CONSERVATION LAW 
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Abstract. Equation equivalents au troisieme ordre pour les schemas de Boltzmann sur 
reseau a vitesse relative a une loi de conservation. Nous etudions la precision formelle des 
schemas de Boltzmann sur reseau a vitesse relative. Ils different des schemas de d’Humieres puisque 
I’ensemble des moments pilotant la relaxation depend d’un champ de vitesse arbitraire fonction du 
temps et de I’espace. Nous traitons de I’asymptotique des schemas a vitesse relative pour une loi de 
conservation : nous etablissons I’equation equivalente du troisieme ordre pour un nombre arbitraire 
de dimensions et vitesses. 


Abstract. We study the formal precision of the relative velocity lattice Boltzmann schemes. They 
differ from the d’Humieres schemes by their relaxation phase: it occurs for a set of moments parametrized 
by a velocity field function of space and time. We deal with the asymptotics of the relative velocity 
schemes for one conservation law: the third order equivalent equation is exposed for an arbitrary 
number of dimensions and velocities. 

Version FRANgAiSE abregee 

Dans cette contribution, nous presentons I’equation equivalente au troisieme ordre des schemas de 
Boltzmann sur reseau a vitesse relative introduits dans |8l |TT] (inspires par le schema cascade |12| 1 
lorsqu’il n’y a qu’une loi de conservation. Ces schemas a vitesse relative different des schemas de 
d’Humieres [5] par leur phase de relaxation : I’ensemble des moments est maintenant parametre par 
un certain champ de vitesse fonction du temps et de I’espace. 

En dimension d G N* d’espace, nous considerons C un reseau cartesien de pas Ax. Le pas de temps At 
est determine par une mise a I’echelle acoustique apres specification d’une vitesse A = Ax/At. Nous 
notons DdQq un schema a q vitesses Vj G 0 ^ j ^ auxquelles sont associees les distributions 
de particules fj, 0 ^ j ^ q—1. 

Le schema consiste a faire evoluer les distributions en deux etapes : la phase de collision puis la phase 
de transport. En notant (Pq, ... ,Pg_i), I’ensemble des polynomes deffnissant les moments, la matrice 
de passage des distributions a ces moments est definie par m ou u{x, t) G est le parametre champ 
de vitesse relatif pour x G t G M. Le vecteur des moments m{u) = {mQ{u),... ,m^_^(u)) est 
alors defini par ([2|) et la phase de relaxation est donnee par ([3]), ou rnj/^{u) est I’equilibre du moment 
rrif^i/u) et s = (sq, .. ., G le vecteur des parametres de relaxation. La phase de transport est 
donnee par ([7]), apres retour aux distributions. Le choix tr = 0 permet de retrouver le schema de 
d’Humieres. 

Nous travaillons avec une seule loi de conservation sur la densite p definie par ([5]) : ainsi 6 q = 0 et 7^ 0, 
1 ^ k ^ q—1. L’objectif est la determination de I’asymptotique du troisieme ordre d’un schema D^Qq 
general. Nous utilisons pour cela la methode des equations equivalentes, un developpement formel 
initialement utilise pour I’etude de schemas aux differences finies [HllS]. Cette methode a ete adaptee 

1 


2 


FRANgOIS DUBOIS, TONY FEVRIER, AND BENJAMIN GRAILLE 


aux schemas de d’Humieres dans mE] et utilisee pour les schemas a vitesse relative a deux lois de 
conservation | 8 ]. Elle s’effectue grace a un developpement de Taylor de I’equation de transport ([7]) a 
petit pas de temps At. Pour cela, les distributions de particules sont supposees etre les restrictions sur 
le reseau £ de fonctions suffisamment regulieres. Les parametres de relaxation sont supposes constants 
et I’equilibre, provenant d’un vecteur de distributions ■ ■ ■, fg^i) independant de u, verifie 

@- 

Afin d’alleger les notations, un indice latin dans les signes de summation decrit implicitement les 
entiers de 0 a q—1, tandis qu’un indice grec decrit les entiers de 1 a d. Nous noterons par ailleurs 
A = At. Sous ces hypotheses, nous communiquons le resultat suivant. 

Proposition L’equation de conservation sur la densite s’ecrit au troisieme ordre grace a I’equation 

m- 

q-l 

La quantite de mouvement a I’equilibre est donnee par derivee particulaire d^ 

j=o 

par dj = 0 ^ j ^ “ 1 ) le default de conservation 0f.{u) par 0j,{u) = d^ 

a j 

0 ^ ^ q — 1, et le parametre de Henon cr^ [H] par <7^ = 1/s^ — 1/2, pour / 0, 0 ^ fe ^ q — 1. 

Ce resultat s’obtient en developpant successivement I’equation de transport ([7]) aux ordres de zero a 
trois. La composition par la matrice des moments donne le developpement equivalent sur les quantites 
macroscopiques. La cle pour monter en ordre est I’utilisation de lemmes de transition m- ces lemmes 
sont les developpements en A a differents ordres des moments non conserves par la collision. 

Dans I’equation de conservation (| 8 ]), nous remarquons que les termes d’ordre inferieur ou egal a un 
sont independants de u grace a I’hypothese (|3]) sur I’equilibre. Nous rappelons que ce constat est aussi 
vrai dans le cas a deux lois de conservation [ 8 ]. Le parametre u apparait dans le terme d’ordre deux : 
il est susceptible d’influer sur le comportement dispersif du schema. Nous notons enfin I’importance 
du defaut de conservation : cette quantite, dependante de I’equilibre, pilote les termes d’ordre un 
(diffusion) et deux (dispersion). 


Introduction 

The analysis of consistency for the lattice Boltzmann schemes uses some Chapman-Enskog type formal 
expansions. This method, based on a multiscale expansion at small Knudsen number |3], has been 
a useful tool for these schemes 13 El- In this contribution, we choose an alternative method: the 
equivalent equations development proposed in dmn] for the study of finite difference schemes and 
adapted to the lattice Boltzmann framework in j 6 ]. The third order equivalent equations of a general 
d dimensions and q velocities——d’Humieres scheme |5] for one and two conservation laws have 
been determined with this method [Tj- It has also been used to build high order schemes [Hilo] , to 
enforce isotropy to an arbitrary order mm and to study some vectorial lattice Boltzmann schemes 
for hyperbolic systems m- 

This method has already been extended to the relative velocity D^Qg scheme for two conserva¬ 

tion laws. We adapt the work to one conservation law: the third order equivalent equation is derived. 
In a first section, we recall the characteristics of the relative velocity D^Qq scheme and exhibit the 
differences with the d’Humieres scheme. The second section introduces the main result of the paper, 
the third order expansion on the conserved variable. The steps leading to this formal development are 
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described in a third section. The influence of the relative velocity is studied: we particularly locate the 
terms involving it. We also check the consistency with the asymptotics obtained for the d’Humieres 
scheme [7]. 


1. The relative velocity D^Qg scheme 

In this section, we recall the derivation for one conservation law of the relative velocity lattice Boltz¬ 
mann schemes introduced in |8] and inspired by the cascaded scheme m- We consider C, a cartesian 
lattice in d dimensions with a typical mesh size Ax. The time step At is determined by the acoustic 
scaling after the specification of the velocity scale A G M by the relation At = Ax/A. For the scheme 
denoted by D^Qg, we introduce V = {vq, ... a set of q velocities of We assume that for 

each node x of C, and each Vj in V, the point x + “VjAt is also a node of the lattice C. The aim of 
the D^Qg scheme is to compute a particle distribution f = {Jq, ■ ■ ■, fg-i) on the lattice C at discrete 
values of time. 

The scheme splits into two phases for each time iteration: first, the relaxation that is non linear and 
local in space, and second, the linear transport for which an exact characteristic method is used. In 
the framework of the relative velocity lattice Boltzmann schemes, the relaxation is written into a basis 
of moments depending on a velocity field u{x,t), a given function of space and time. In the following, 
we note the velocity field u because it appears only during the collision that is local in space. The 
matrix of moments M{u), which is supposed to be invertible, is defined by 

^(^)kj = Pki^j ( 1 ) 
where {Pq, ..., P^_^) are some polynomials of M[A^,..., A^]. We choose the first polynomials as 
follows: ig = 1, = Xj ^, 1 ^ k ^ d. The moments m{u) = (mQ(xt),... are then given by 

m{u) = M{u)f, (2) 

so that m^iu) is the moment in the frame moving at the velocity u. The relaxation phase is 
diagonal in the shifted moments basis 

m^(u) = 0 ^ k ^ q-1, (3) 

where rn^{u) is the k^^ moment at equilibrium and Sp., the relaxation parameter associated with the 
k^^ moment for 0 ^ k ^ Q—^- The vector of the distribution functions at the equilibrium f'^'^ is chosen 
independent of the velocity field u so that we have 

m^^{u) = M{u)f^^ = Miu)M{0)-^m^^{0). (4) 

We choose to work in the case of one conservation law on a scalar variable named here “density” 

defined by 

P = E4- w 

3 

Consequently, the associated relaxation parameter Sg is null. The particle distributions are then 
computed by 

f* = M~^ {u)m* (u). (6) 

The transport phase reads 

fj{x,t + At) = f*{x - v-At,t), 0 ^ j ^ q-1. 

Let us note that the d’Humieres scheme is enclosed in this framework by taking u = 0. 


( 7 ) 
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2. Third order equivalent equation 

A Taylor expansion method |7j is performed to derive the third order equivalent equation of the 
relative velocity D^Qq scheme for one conservation law, with any general equilibrium. We explicit a 
formal development of the scheme when At and Ax go to zero with A = Ax/At and the relaxation 
parameters s fixed. We obtain in this way a set of partial differential equations that are consistent 
with the scheme at different orders. 

In the following, we do not adapt the Boltzmann scheme to a particular partial differential equation— 
at first or second order—but we put in evidence the various operators hidden inside the algorithm. The 
reasoning consists in a formal development of the transport phase ([7]) at small At, assuming that all 
particle distributions are the restrictions on the discretized space of a sufficiently regular distribution 
function. The Taylor formula can thus be used as much as wanted. 

The result we want to emphasize is given by the proposition 12.11 this is the third order equivalent 
equation for the relative velocity scheme with one conservation law. We particularly prove that 

the second order equations are the same as the ones of the d’Humieres scheme [3: this indepen¬ 
dence results from our choice for the equilibrium Q. The relation (|3]) means that the approximated 
physics do not depend on the velocity field u up to the second order. The dependence only appears 
at the third order. In the following, we note A = At. Unless otherwise specified, a sum over a greek 
parameter goes from 1 to d and over a latin parameter from 0 to — 1. 

Proposition 2.1. The third order conservation equation on the density reads: 

^tP + ^ u/ 9^(d^ (M(S)-/0;(S)) 

/3 p 

+ ^ E 'iT"D + ^ E (8) 

q-l 

where the momentum equilibrium is given by particular derivative by 

j=0 

dj = 4 + 0 ^ j ^ q — 1, the conservation default 0^{u) by 9^ {u) = . d^ fp, 

0 ^ A; ^ g — 4, and the Henon’s parameter |14) by a^ = \ls^ — 1/2, for Sj^ ^ 0, 0 ^ k ^ q — 1. 

We observe that each term of order I in A contains only space derivatives of order I + 1. 

3. Principal ideas for the proof of the main result 

The detailed calculations of the proof of the proposition 12.11 are available in m To obtain the 
equivalent equation up to any order p G N*, a multidimensional Taylor formula relative to A is applied 
to the transport phase ([7]). This expansion is then written in the moments basis after multiplying by 
M{u)kj 

/\l _ ('_Ail“l 

E E ^-W(",)"v(S)t,A/; + o(Af+i), o^k^q-1, (9) 
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d d d 

where a = (ai, ... ,ad) € N'^, : 9^ = da^ ■ ■ ■ a! = JJoq,! and |o| = ^ a«. 

a.=l a.=l Q;=l 

This relation develops the different moments according to the choice of the moment index k. Taking 
/c = 0 yields to the conservation of the density, and A: > 1 to the expansions of the non conserved 
moments. 

To obtain the third order equivalent equation on the density, we proceed order by order. The first 
step thus consists in deriving the zeroth order corresponding to p = 0 in ([9]). This gives the following 
proposition: its proof can be found in [6j. 

Lemma 3.1 (The particles stay close to the equilibrium). The particle distributions can be 
expanded as 

/^. = /;'! +0(A), /; = //‘i + 0(A), Q^j^q-l. 

These identities are assumed to be formally derivable as much as wanted: this point is used to obtain 
the first order equivalent equation on p taking p = 1 and A: = 0 in ([9]) . This equation is the truncation 
of ([8]) at the first order. 

The key to gain an order is to develop relatively to A the non conserved moments. To get the second 
order equivalent equation on p, these moments must be expanded at the second order. We present 
these expansions in the third order case, the second one being obtained by truncation. 

Lemma 3.2 (Third order transition lemma). The non conserved moments before and after the colli¬ 
sion read, 


m^{u)= m^‘^(w)-A(^^ 

+ A,<t) 

+ 0(A3), 

1 ^ k ^ q - 1, 

(10) 

m*(w) = m®'‘(S) + A(^i 

b 

<r 

b 

1 

+ 0(A3), 

1 ^ k ^ q - 1, 

(11) 


fk{u,A,a) = 0fc(w)-A ai M{u)^. dl (M(w)^./0;(w)). 


Two relations are combined to obtain these identities: the expansion ([9]) at the desired order and the 
relaxation phase Q . They give immediately the second order truncation of the relations (|10llll) . This 
second order version is the keypoint to get the complete third order (|10llll) . 

To get the second order equivalent equation on the density, we choose p = 2 and A = 0 in ([9]). The 
calculation necessitates the second order expansion of the post collision momentum that is given by 
IfTT]) taking 1 ^ k ^ d. 

To get the complete third order equivalent equation on p, the same reasoning is used an order further. 
The required third order development of the post collision momentum is obtained using (ttH) with 
1 ^ k ^ d. The second order truncation of (|8]) is then useful to characterize this development up to 
the third order, that leads to the proposition 12.11 

Let us now recover the equation on the density for the d’Humieres scheme |7] corresponding to S = 0. 
We focus on the second order terms, the first order being independent of u. The quantity 

0 
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is already present in [7]. The second term, independent of the relaxation parameters, reads 

^ vfv]dl^{dl fp) = dP{Ap{0)di{0)). 

Pa,I 

There is an only term depending on cr in (jH]) . Splitting the time and space derivatives of the particular 
derivative d^, it reads 

^ vf dpdl (M(O)-^0;(O)) = Y + Y. 

where Ap{0) is the momentum velocity tensor defined by A^,^(0) = E vP^M{0)-\ 1 ^ /3,7 < d, 
0 ^ ^ q — 1. We have thus recovered all the terms exhibited in [7]. j 
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